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Summary 

System-identification methods compose a mathematical 
model, or series of models, from measurements of inputs 
and outputs of dynamic systems. The extracted models 
allow the characterization of the response of the overall 
aircraft or component subsystem behavior, such as actua- 
tors and on-board signal processing algorithms. This 
paper discusses the use of frequency-domain system- 
identification methods for the development and integra- 
tion of aircraft flight-control systems. The extraction and 
analysis of models of varying complexity from nonpara- 
metric frequency-responses to transfer-functions and 
high-order state-space representations is illustrated using 
the Comprehensive Identification from ErEquency 
Responses (CIFER®) system-identification facility. 
Results are presented for test data of numerous flight and 
simulation programs at the Ames Research Center includ- 
ing rotorcraft, fixed-wing aircraft, advanced short takeoff 
and vertical landing (ASTOVL), vertical/short takeoff and 
landing (V/STOL), tiltrotor aircraft, and rotor experiments 
in the wind tunnel. Excellent system characterization and 
dynamic response prediction is achieved for this wide 
class of systems. Examples illustrate the role of system- 
identification technology in providing an integrated flow 
of dynamic response data around the entire life-cycle of 
aircraft development from initial specifications, through 
simulation and bench testing, and into flight-test 
optimization. 

1. Introduction 

System identification is a procedure for accurately charac- 
terizing the dynamic response behavior of a complete air- 
craft, subsystem, or individual component from measured 
data. This key technology for modem fly-by-wire flight- 
control system development and integration provides a 
unified flow of information regarding system performance 
around the entire life cycle from specification and design 
through development and flight test (fig. 1). A similar 
“roadmap” for application of system-identification meth- 
ods to rotorcraft development was previously proposed by 


Schrage in a comprehensive report dedicated to this topic 
(Anon. 1991). An excellent historical summary and 
overview of system identification is given by Hamel 
(1995). System identification has been widely utilized in 
recent aircraft programs including many of those 
described in the present volume. Common applications for 
flight-control system development include: definition of 
system requirements, specification and analysis of han- 
dling qualities, evaluation of proposed control-law con- 
cepts, validation and improvement of complex simulation 
models, validation of subsystem components and devel- 
opment facilities, and flight-test optimization of control 
laws. 

Frequency-domain identification approaches are espe- 
cially well suited to the development and validation of 
flight-control systems. Feedback stability and noise 
amplification properties are determined from the broken- 
loop frequency response, and characterized by metrics 
such as crossover frequency, and associated gain and 
phase margins. Command tracking performance is deter- 
mined from the closed-loop frequency-response, and 
characterized metrics such as bandwidth and time-delay, 
and equivalent system eigenvalues. Frequency-domain 
identification approaches allow the direct and rapid 
(including real-time) identification of these frequency 



Figure 1. Road map for fly-by-wire flight-control system 
development and integration. 






responses and metrics, without the need to first identify a 
parametric (state-space) model structure such as is 
required in applying time-domain methods. Careful 
tracking of the broken-loop and end-to-end closed-loop 
frequency-response behavior from the preliminary design 
studies through detailed design and simulation and into 
flight test provides an important “paper trail” for docu- 
menting system performance and solving problems that 
may appear in the later phases of development. 

The availability of comprehensive and reliable computa- 
tional tools has substantially enhanced the acceptability of 
frequency-domain techniques in the flight-control and 
flight-test community. Benefits from applying these tech- 
niques include the reduction of flight-test time required 
for control system optimization and handling-qualities 
evaluation, especially for complex control-law architec- 
tures, and improvements in the final system performance. 
Frequency-domain methods offer a transparent under- 
standing of component and end-to-end response character- 
istics that can be critical in solving system integration 
problems. 

This paper reviews frequency-domain system- 
identification methods for development and integration of 
aircraft flight-control systems. These methods were 
developed under a long-term research activity at the Ames 
Research Center by the Army Aeroflightdynamics Direc- 


torate (AFDD), the NASA and Sterling Software. Many 
of the flight applications have been to rotorcraft, which 
pose an especially difficult challenge to system identifica- 
tion (Tischler 1990). The dynamics of these aircraft are 
highly coupled, and unstable. Additionally, the rotorcraft 
dynamics include lightly-damped fuselage and rotor 
modes. Vibration and low excitation signal content, espe- 
cially near hover, results in typically low signal-to-noise 
ratios. Experience in developing and applying system- 
identification methods to the rotorcraft problem has pro- 
duced a set of tools that has proven highly reliable for the 
broad scope of applications reviewed in this paper. The 
first section presents a summary of the frequency-domain 
approach and the Comprehensive Identification from 
FrEquency Responses (CIFER®) comprehensive analysis 
facility. The remainder of the paper is organized into five 
sections following the flight-control development 
flowchart of figure 1 from specifications and design 
through flight-test optimization. Each section illustrates 
important techniques with examples based on fixed and 
rotary-wing projects at the Ames Research Center. 

2. Overview of AFDD/NASA System- 
Identification Techniques 

The AFDD/NASA frequency-domain system identifica- 
tion procedure is shown in figure 2, and in reviewed in 



Figure 2. Frequency-domain system identification procedure. 
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this section. Details of the procedure are found in Tischler 
and Cauffman (1992). System-identification methods and 
requirements for specific application to high-bandwidth 
rotorcraft flight-control system design are given by 
Tischler (1990). 

Aircraft or subsystem component dynamics are excited by 
a pilot-generated or computer-generated frequency-sweep 
input. The dynamic responses are generally measured by 
dedicated sensors, and the data are either recorded on- 
board or telemetered to the ground for processing. 

Kalman filtering techniques (or simple numerical integra- 
tion) are used to check data compatibility and eliminate 
spurious instrumentation system biases, scale factors, and 
drop-outs. Here, unmeasured signals may be estimated 
from the available measured states. 

The foundation of the AFDD/NAS A approach is the high- 
quality extraction of a complete multi-input/multi-output 
(MIMO) frequency-response database. These responses 
fully characterize the linearized coupled characteristics of 
the system without a priori assumptions. Advanced multi- 
variable spectral analysis using the Chirp-Z transform and 
composite optimal window techniques have been devel- 
oped and exercised over many flight applications. These 
methods provide significant improvement in identification 
quality relative to standard fast Fourier transform (FFT) 
methods. The frequency-response database directly 
supports important flight-control system applications 
including: handling-qualities analysis and specification 
compliance testing, simulation validation, and servo-loop 
stability analysis. 

Transfer-function fitting is a rapid procedure for extract- 
ing simple single-input/single-output parametric models 
of specific frequency-responses pairs. These transfer- 
function models define the lower-order equivalent sys- 
tems (LOES) of the fixed wing handling-qualities specifi- 
cations (MIL-STD-1797) and directly support root-locus 
techniques for flight-control system design. 

Accurate MIMO state-space models are often needed to 
support multivariable control-law design, simulation 
model validation and improvement, and validation of 
aerodynamic theory or wind tunnel results. Here, sophisti- 
cated nonlinear search algorithms are used to extract a 
general state-space model that matches the complete 
MIMO input/output frequency-response database. A sig- 
nificant advantage of identifying parametric models from 
frequency responses is the capability to individually 
define the appropriate frequency range for each response 
pair based on the associated coherence function — a valu- 
able accuracy and linearity metric. The coherence func- 
tion is also useful for automatically selecting error 
weighting in the cost function independent of the model 
structure. A methodical and integrated model structure 


determination procedure simplifies the model to a mini- 
mum set of reliable parameters that accurately character- 
izes the MIMO frequency-response database. Finally, the 
identified state-space model is validated by comparing 
predicted time responses with the actual flight responses 
for test inputs not used in the identification procedure. 

The frequency-domain system-identification procedure is 
incorporated in a sophisticated interactive computational 
facility known as CIFER®-Comprehensive Identification 
from ErEquency Responses. Integrated data-basing and 
extensive user-oriented utilities are distinctive features of 
CIFER® for organizing and analyzing the large amounts 
of data which are generated for flight-test identification 
projects. A screen-driven interface is tied to the database 
for rapid user interaction. Previous program set-ups and 
analysis results are retrieved by simply referencing case 
names. Then, changes can be easily made by moving the 
cursor around on the user screens and modifying the 
default or previously saved program parameter values. 

The changes are then updated in the database with a single 
key stroke. Utilities are available for quick inspection, 
searching, plotting, or tabulated output of the contents of 
the database. Extensive analysis modules within CIFER® 
support: 1) rapid identification of transfer- function mod- 
els; 2) signal spectral analysis; 3) handling qualities and 
classical servo-loop analysis; and 4) time and frequency- 
domain comparisons of identification and simulation 
model predictions with flight data. Aircraft applications of 
CIFER® have included the full life-cycle of flight-control 
system development depicted in figure 1. The Deutsche 
Forschungsanstalt fur Luft- und Raumfahrt (DLR) Insti- 
tute for Flight Mechanics (Braunschweig, Germany) has 
also developed and widely applied excellent methods for 
frequency-domain system identification. Applications to 
flight-mechanics and flight-control studies at the DLR 
include rotorcraft, transport aircraft, and high- 
performance aircraft (Kaletka and von Grunhagen 1989, 
Kaletka and Fu 1993). 

3. Design Specifications and Specification 
Acceptance Testing 

Formulating design specifications is the starting point for 
flight-control system development, while validating the 
achievement of these design goals is the concluding step 
in the process (fig. 1). Dynamic models of expected sys- 
tem behavior are determined in the design process using 
system identification and are tracked and updated 
throughout the aircraft development and flight testing. 
This documentation provides an important “paper trail” 
that minimizes flight-control development time and 
reduces the need for costly flight-test tuning. This section 
presents system-identification methods for defining and 
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verifying design specifications. Flight test examples illus- 
trate the analysis of handling qualities and servo-loop sta- 
bility characteristics. 

Early handling-qualities specifications for fixed-wing air- 
craft (MIL-F-8785A, Anon. 1954), and for rotary- wing 
aircraft (MIL-H-8501A, Anon. 1961) were based on sim- 
ple dynamic modeling concepts and time-domain metrics. 
These specifications were suitable to aircraft in which 
stability augmentation systems (SAS) did not significantly 
alter the character of the (classical) bare-airframe flight- 
mechanics responses. Compliance testing techniques 
depended on standard step and doublet inputs long used in 
the flight-test community, with little requirement for 
sophisticated post-flight-data processing. 

Modem fly-by-wire aircraft employ high-bandwidth digi- 
tal flight-control systems to achieve greatly increased 
agility and disturbance rejection across a significantly 
widened operational flight envelope as compared with the 
older generation of aircraft. The flight control includes 
complex feedback and feedforward shaping and advanced 
control moment devices that profoundly alter the bare- 
airframe characteristics and invalidate the classical stabil- 
ity and control modeling concepts and testing methods. 
For example, modem combat aircraft achieve independent 
pitch pointing and flightpath control with direct lift 
devices and vectored thrust, rather than the coupled 
attitude-path response to elevator for conventional air- 
craft. This capability greatly enhances weapon pointing 
and air-to-air combat maneuvering. Another common fea- 
ture of advanced aircraft is side-stick controllers which 
reduce weight, space, and cockpit complexity compared 
to standard center-sticks. Classical static stick-stability 
testing is an invalid method for determining speed stabil- 
ity since the side-sticks possess automatic trimming at 
neutral stick position and feedback loops provide the 
required stability independent of the trim gradient. 

A new concern that arises for modem fly-by-wire aircraft 
is the potential for the accumulation of effective time 
delays due to digital flight-control computations, flight- 
control system filters, and fly-by-wire actuators. Actuator 
rate-limiting can also contribute large equivalent time 
delays in modem aircraft (Buchholz et al. 1995). Exces- 
sive delays have been repeatedly cited as a key cause for 
handling-quality problems and stability-loop margin 
degradation in modem aircraft, yet equivalent time delay 
can not be reliably measured using the standard testing 
techniques. Clearly the dynamics modeling concepts, 
specifications, and testing techniques must be appropriate 
to the unique characteristics of modern highly-augmented 
aircraft. 

System identification provides an accurate, rapid, and 
reliable approach for defining design specifications and 


for validating aircraft flight performance for highly- 
augmented flight-control systems. The modem U.S. fixed- 
wing specification (MIL-STD-1797; Anon. 1987) and 
rotorcraft specification (ADS-33C; Anon. 1989) are based 
on extensive frequency-domain system-identification 
analyses of flight-test and simulation responses. Numer- 
ous examples from these and comparable European 
handling-qualities specifications are presented in this vol- 
ume and in the references of this paper. Two common 
handling-quality specifications are the bandwidth/phase- 
delay criteria and the LOES criteria. The former is 
checked directly from frequency-response identification, 
and the latter is checked from a transfer function fit of the 
frequency-response result. An illustration of flight-test 
and handling-qualities analyses based on these specifica- 
tions is now presented. 

The Advanced Digital Optical Control System (ADOCS) 
demonstrator (fig. 3), developed by Boeing’s Helicopters 
Division under contract to the U.S. Army, was a UH-60A 
helicopter highly modified with redundant processors, 
instrumentation, and side-stick controllers (Glusman et al. 
1987). The overall program objective of the ADOCS was 
to provide the technology base for the engineering devel- 
opment of an advanced battlefield-compatible flight-con- 
trol system that: 1) enhanced aircraft mission capability; 

2) improved handling qualities; and 3) decreased pilot 
workload. System identification flight tests and analyses 
using CIFER® were conducted to document the response 
characteristics and to compare handling-qualities charac- 
teristics with the (proposed at that time) ADS-33 design 
specifications (Hoh et al. 1988). Aircraft excitation was 
achieved via piloted frequency sweeps using the side-stick 
controller as shown in figure 4. Real-time telemetry of 
pilot inputs and aircraft responses ensured that pre- 
established aircraft flight limits were not exceeded. 

The ADOCS frequency response for pitch response due to 
longitudinal input is shown in figure 5, along with the 



Figure 3. Advanced digital optical control system (ADOCS) 
demonstrator. 
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Time (sec) 

Figure 4. Longitudinal side-stick frequency-sweep in hover; 
(a) pitch input, (b) pitch rate. 



Frequency (rad/sec) 


Figure 5. Identification of ADOCS pitch-rate response in 
hover. 


determination of bandwidth and phase delay as required 
by the ADS-33 specification. The value of the coherence 
function (fig. 5) is consistently above 0.8 for the fre- 
quency range of 0.2-8 rad/sec indicating excellent identi- 
fication. At higher frequencies, the coherence drops, 
which reflects the intentionally reduced piloted inputs. 

The pitch bandwidth and phase delay values obtained 
from the hover identification results are shown on the 
ADS-33 specification boundary in figure 6. Level I (satis- 
factory) handling qualities for moderate pilot-gain tasks 
such as the helicopter “boh-un” are predicted, which is 
consistent with the Cooper-Harper pilot rating displayed 
next to the data symbol. A good correlation of pilot rating 
and the predicted handling qualities result was also indi- 
cated in pitch for the 80 kt flight condition. 

Transfer-function (LOES) models of the key on-axes 
responses were generated from single-input/single-output 
fits of the identified ADOCS frequency responses. The 
LOES pitch response in hover is: 

9 - 0,876 (s + 0,229) e~°' 238s 

SLON " s [0.539, 1.82] 

The response is dominated by a well-damped second- 
order mode with a frequency of 1.8 rad/sec. The LOES 
handling-qualities specification boundaries of figure 7 
have been established based on system-identification 
analyses of an extensive flight-test and simulation data 
base (Hoh and Ashkenas 1979). The ADOCS characteris- 
tics are seen to lie well within the Level I region. A sec- 
ond important characteristic of the LOES pitch response is 
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“BWq ( rac ^ sec ) 

Figure 6. Handling-qualities correlation of ADOCS small 
amplitude pitch response in hover and 80 kt; average pilot 
ratings are shown next to the data. 
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Equivalent time delay, x 0 (sec) 


Figure 7. Correlation ofADOCS LOES pitch response with 
handling-qualities data. 

the relatively large equivalent time delay, x = 238 msec. 
Handling-qualities experience indicates that the equiva- 
lent time delay should not greatly exceed % = 120 msec, 
thereby suggesting ADOCS handling-qualities degrada- 
tion for “high-gain” tasks. Comparable levels of equiva- 
lent time delay in the roll axis were considered to be a key 
contributor to pilot-induced oscillations (PIO) for “high- 
gain” piloting tasks such as slope landing (Tischler et al. 
1991). Additional examples of handling-qualities analyses 
using lower-order equivalent-system modeling are pre- 
sented in this volume. 

The stability characteristics of aircraft rigid-body and 
structural dynamics may be greatly affected by the feed- 
back loops of the flight-control system. Feedback may 
degrade the flutter margin stability at the same time it 
improves the rigid-body stability and handling qualities. 
Military specification 9490D (Anon. 1975) defines mini- 
mum levels of control system gain and phase margin as 
determined from a broken servo-loop frequency-response 
analysis. The specifications are given as a function of fre- 


quency range, with larger gain margins required for the 
structural elastic modes (table 1). Figure 8 shows the bro- 
ken servo-loop frequency response of a large single-rotor 
helicopter as obtained by computer-generated frequency- 
sweep flight-test procedures. The rigid-body response 
crossover frequency is 2.0 rad/sec with an associated 
phase margin of 28 deg. This margin is slightly below the 
recommended specification value. The gain margin is 
checked at each crossing of the 180 deg (±n360) phase 
line as shown in the figure. The critical margin is the min- 
imum value (GM 5 ), which is 15 dB at a frequency of 
23.5 rad/sec. This frequency corresponds to the first verti- 
cal bending mode for the tail boom of this aircraft. 
Reference to table 1 indicates this gain margin to be well 
within accepted design specifications. The coherence 
function of figure 8 shows excellent identification accu- 
racy for this flight test across the broad frequency range of 
interest (1-30 rad/sec). Sharp drops (“holes”) or peaks in 
the coherence function reflect structural anti-nodes and 
nodes respectively. Examples of fixed-wing programs 


Table 1. MIL-F-9490D gain and phase margin requirements (dB, deg) 
from Caldwell (1994) 


Airspeed 
Mode Freq. Hz 

Below 

^omhn 

Vomin 

to 

v 

v o max 

At 

Limit Speed 

V L 

At 

1.15 *V L 

f M < 0.06 

GM = 6dB 
No Phase 
Reqt. 
below 
V 

v omin 

GM = ± 4.5 
PM = ± 30 

GM = ± 3.0 
PM = ± 20 

GM = 0.0 
PM = 0.0 
Stable 
at 

Nominal 

Phase 

and 

Gain 

0.06 £ f M < tStASE 

GM = ± 6.0 
PM = ± 45 

GM = ± 4.5 
PM = ±30 

1st ASE < 

GM = ± 8.0 
PM =±60 

GM = ± 6.0 I 
PM = ±45 


Stability-loop: AVCS-output/servo-command input 



Frequency (rad/sec) 


Figure 8. Rotorcraft sen/o-elastic stability margin testing. 
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using frequency-domain system identification for elastic 
mode stability -margin evaluation include the X-29 
(Clarke et al. 1994) and EAP (Caldwell 1994) aircraft. 

4. Design 

The design process establishes control system loop archi- 
tecture and associated control-law parameters that achieve 
desired handling-qualities and servo-loop stability specifi- 
cations. During the conceptual design phase, system- 
identification procedures are applied to simple linear- 
design models to establish a baseline description of the 
proposed control system approaches and an initial check 
of specification compliance. At the detailed design stage, 
system-identification methods can extract highly-accurate 
linear-control system design models from complex simu- 
lation models or wind-tunnel data. These applications of 
system identification in the design process are illustrated 
in this section. 

Conceptual control -system design studies are commonly 
based on simple stability and control derivative descrip- 
tions and transfer function of the airframe dynamics as 
obtained from first-principles aerodynamic theory. 
Control-law architectures are conceived and the initial 
system is modeled in a computer-aided design (CAD) 
environment such as MATLAB® (1992). System identifi- 
cation provides LOES models which are very useful in 
characterizing the end-to-end system dynamics and 
delays, and for an initial check against the design specifi- 
cations. Flight-control system design parameters are then 
adjusted until the identified LOES characteristics satisfy 
the design requirements. In the F- 15 S/MTD demonstrator 
project, a numerical optimization design tool was devel- 
oped to automatically adjust control-law parameters to 
meet LOES specifications (Moorhouse and Citurs 1994). 

Detailed flight-control design efforts are based on very 
complex high-order and nonlinear-simulation models. 
Force and moment descriptions are developed for each of 
the aircraft elements such as the wings, propulsion sys- 
tem, and flight-control systems based on wind tunnel 
look-up tables, component bench-test data, and analytical 
theory. The simulation of multiple rigid-body systems, or 
flexible bodies involves sets of dynamic equations of 
motion linked by constraint conditions. In many simula- 
tions these sets of equations are numerically integrated in 
serial form to reduce the complexity of deriving a fully- 
coupled multibody simulation. The distributed or serial 
nature of these complex simulations thus may preclude 
the extraction of an integrated high-order linear model of 
the fully-coupled system as is needed for accurate control- 
design studies. 


Even when the simulation architecture allows for the 
direct extraction of higher-order linear models using clas- 
sical numerical perturbation methods, the assumption of 
independent perturbations results in incorrect phasing of 
the state variables within the multidimensional look-up 
tables. For example, the look-up table for aerodynamic 
pitching moment may depend both on angle of attack and 
pitch rate, so = f(a). Thus, the correct determination 
of phugoid dynamics depends on maintaining representa- 
tive phasing of q and a within the linearization process. 
Selection of perturbation size can also strongly influence 
the linearization results. These effects can significantly 
degrade the predictive accuracy of the extracted linear 
model. Much more accurate linear models are obtained by 
simulating piloted frequency-sweep inputs and extracting 
state-space models using system identification just as if 
from flight-test data. 

Engelland extracted accurate stability and control 
derivative models of a conceptual A/STOVL air- 
craft from a complex nonlinear off-line simulation 
(Engelland et al. 1990) to support control-system 
design studies. The excitation input consisted of 
computer-generated frequency sweeps and white 
noise. In a procedure described by Ballin and Dalang- 
Secretan (1991), artificial feedback control loops 
were included to keep the aircraft flight condition 
near the reference trim point during the inputs. Start- 
ing from the perturbation derivative results CIFER® 
was used to identify a more accurate 6 DOF bare- 
airframe model. The perturbation and CIFER® deriva- 
tives are compared in table 2. Longitudinal-frequency 
responses of the two linear models are compared with 
the complete simulation responses in figure 9 for a 
flight condition of 120 knots. The linear model 
obtained using system identification is seen to be 
much more accurate than the numerical perturbation 
model for the high-frequency (3.0-20 rad/sec) pitch-rate 
response q/8e, and for the low-frequency (0. 1-1.0 rad/sec) 
longitudinal-acceleration response a x /6e. The models 
are essentially identical in the mid- frequency range. 

A time domain comparison of the two-linear models 
with the nonlinear-simulation response is shown in 
figure 10 for a small (1 deg) pitch-doublet input. The 
system-identification model is seen to track the non- 
linear behavior much more closely than the numerical 
perturbation model. The improvements are most notice- 
able for the long-term response (low-frequency) 
behavior, which is consistent with the frequency- 
response comparison of figure 9. 
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Table 2. Comparison of ASTOVL perturbation derivatives and 
CIFER® results 


Derivative 

Perturb. Value 

CIFER Value 

C.R. (%) 

Insens. (%) 


-0.03471 

-0.03602 

-5.662 

2.289 

X w 

0.03958 

0.02852 

6.910 

2.840 

x vv 

6.764E-04 

6.764E-04t 

— 

— 

x <? 

0.2451 

0.2451 f 

— 

— 

X PCD 

-7.690E-03 

-8.303E-03 

-7.504 

3.584 

X PM 

0.02270 

0.02229 

3.731 

1.835 


-0.5150 

-0.5586 

-2.353 

1.005 


-0.04596 

-0.03312 

-13.62 

4.579 


-0.3704 

-0.2817 

-4.386 

1.377 

Z w 

-0.01023 

-0.01 023t 

— 

— 

z q 

-3.754 

-3.754t 

— 

— 

Z PCD 

0.1389 

0.1551 

5.571 

2.698 

ZpM 

-0.3800 

-0.3305 

-2.254 

1.016 

Z©M 

-0.01724 

-0.03055 

-4.646 

2.242 

M u 

1.661E-04 

-1.059EO3 

-6.016 

1.745 

M w 

1.222E-03 

3.71 5E03 

5.263 

1.283 

Mvv 

-1.286E-03 

-1.286E03t 

— 

— 

M q 

-0.4971 1 

-0.6852 

-5.561 ! 

1.873 

MpcD 

0.02494 

0.02818 

2.517 

0.9822 

m PLA 

4.993E-04 

4.993E04f 

— 

— 

Mq n 

2.502E-04 

4.953E-04 

10.16 

4.257 


t Perturbation value used. 



Nonlinear simulation 
Perturbation linear model 
CIFER linear model 



Time (sec) 


Time (sec) 


Figure 10. Time-domain comparison of perturbation and 
identification models with nonlinear simulation of ASTOVL 
aircraft : 




frequencies where the dynamic responses are larger, and 
correct phasing of the representative motion variables for 
entry into the multidimensional look-up tables is impor- 
tant. A cursory time-domain comparison of the numerical 
perturbation results with the nonlinear-simulation 
response would suggest the presence of strong nonlineari- 
ties in the A/STOVL aircraft dynamics. However, the 
very close agreement of the system-identification model 
with the nonlinear simulation shows that the method of 
linear model extraction is much more important in this 
case than the nonlinear characteristics of the simulation. 


Figure 9. Frequency-response comparison of perturbation 
and identification models with nonlinear ASTOVL 
simulation. 

These results show that system identification provides an 
A/STOVL linear model that will be much more accurate 
than models extracted using numerical perturbation meth- 
ods. The improvement obtained by “flying” the 
frequency-sweep input is especially apparent at low 


Success in achieving maximum control-system perfor- 
mance and robustness in flight depends heavily on the 
predictive accuracy of the linear-design models. The sys- 
tem identification approach provides highly-accurate 
design models for design at specific flight conditions, but 
it is clearly more time intensive than the simple numerical 
perturbation method. This is not a practical approach for 
checking control system behavior at the tens or hundreds 
of off-nominal conditions. 
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5. Simulation 

The detailed implementation of the control-system design 
is evaluated in comprehensive real-time piloted simulation 
trials. System-identification techniques are exercised to 
validate the real-time math model implementation of non- 
linear digital control laws. Also, these techniques are used 
to document simulator motion and visual systems. Once 
flight-test data are available, system identification pro- 
vides an important tool for validating and updating the 
simulation math models. This section illustrates system- 
identification techniques for validating simulation math 
models and simulator validation using XV- 15 tilt-rotor 
(airplane mode) and UH-60A helicopter results. A com- 
panion frequency-domain format is proposed for specify- 
ing simulation model fidelity for the on-axis responses. 
Finally, an analysis based on an A/STOVL piloted simula- 
tion study illustrates the use of system identification for 
determining actuator authority requirements. 

Direct frequency-response comparison of the end-to-end 
performance of the complex simulation model with the 
conceptual design models and specifications constitutes 
an important “dynamic check” which often exposes unex- 
pected processing delays such as in the numerical integra- 
tion procedures, or errors in the digital (Z-plane) 
implementation of control laws. This technique is also 
useful in exposing degradation in control system perfor- 
mance due to high-order structural or other hardware 
dynamics modeled in the advanced design simulation 
model that may not have been taken into account in the 
conceptual studies. 

Simulator visual and motion systems should track the 
math-model response as accurately as possible to ensure 
that the pilot’s cueing environment is correct and that the 
handling-qualities evaluation obtained in the simulator 
reflect what may be expected in flight. Nonlinear com- 
pensation algorithms have been developed by McFarland 
(1988) that offset visual system delays, thereby minimiz- 
ing the mismatch between the simulator visual system 
response and the math model. In work reported by 
Atencio (1993) and illustrated in figure 11, there is nearly 
perfect agreement of the DIG-1 visual system image (with 
McFarland compensation) and the UH-60A helicopter 
simulation math model. Math-model commands to the 
simulator motion drive are attenuated using wash-out 
logic. The wash-out parameters are selected to preserve 
the dynamic behavior in the frequency range of most con- 
cern to the pilot (e.g., 1-10 rad/sec for pitch and roll 
tasks), while accommodating the restricted motion envi- 
ronment of the simulator. Figure 12 from Atencio (1993) 
compares the (washed-out) cab roll-motion to stick input 
with the math-model response for the UH-60A simulation 
in the Ames Vertical Motion Simulator (VMS). In the 



Frequency (rad/sec) 

Figure 11. Validation of DIG-1 visual-system response. 



Frequency (rad/sec) 

Figure 12. Documentation of VMS motion follow-up for 
UH-60A roll response in hover. 

1-10 rad/sec frequency range, the simulator motion drive 
response follows the math model, although the motion is 
less than one-to-one as seen by the vertical shift in the 
magnitude curves. The motion drive wash-out logic is 
designed to minimize phase distortions in this frequency 
range as can be seen in the figure. At low frequencies the 
large motion is washed-out, and considerable errors are 
encountered in the magnitude and phase response as 
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expected. At high frequency, the motion drive is unable to 
follow rapid commands to the aircraft model, resulting in 
larger phase lags of the motion follow-up as seen in the 
figure. 

Once flight-test data are available, system-identification 
tools are exercised to validate and update the simulation 
math models. The direct comparison of frequency- 
response behavior provides a clear picture of model 
fidelity as a function of frequency. This is critical for val- 
idating piloted simulations since the requirements on pilot 
cueing accuracy are also frequency dependent. Tne sepa- 
rate display of the magnitude and phase responses allows 
the sources of simulation discrepancies to be more easily 
determined. For example, an excessive time delay (t) in 
the simulation math model or hardware causes a linear 
phase shift with frequency (<|> = -TCO). Scaling errors in the 
simulation model appear as a clear vertical shift (in dB) in 
the magnitude curve. These effects are all combined in the 
time-domain and therefore are not easily discernible in the 
traditional time-response comparison methods for valida- 
tion. Further, the procedure of overlaying time histories is 
often not very accurate since the flight responses rarely 
begin in a trim quiescent condition. 

Tischler (1987) conducted an extensive flight-test pro- 
gram and simulation math-model validation study on the 
XV- 15 tilt-rotor aircraft shown in figure 13. This tilt-rotor 
math model is based on comprehensive look-up tables of 
full-scale wind-tunnel test data, and detailed theoretical 
models of the rotor-system behavior and rotor-on-airframe 
aerodynamic interference effects. Figure 14 compares the 
flight and simulation roll responses for a flight condition 
of 170 kts. Excellent dynamic response fidelity is seen in 
the close match of the simulation prediction and the 
measured flight response. Figure 15 replots these results 
in terms of magnitude and phase errors as a function of 
frequency. Here 0 dB magnitude and 0 deg phase indicate 
perfect tracking of the flight and simulation results. Also 
shown in the figure are math model mismatch boundaries 
proposed herein for the highest fidelity training simula- 
tions (FAA Level D). These boundaries correspond to the 
LOES mismatch criteria from the fixed-wing handling- 
qualities criteria (Hoh et al. 1982). The XV-15 simulation 
math model complies with the proposed Level D (high- 
fidelity) criteria. This result is consistent with the very 
favorable pilot comparison of simulator and flight behav- 
ior (Churchill and Dugan 1982). The same approach has 
also been proposed and applied by the DLR for evaluating 
the fidelity of in-flight simulation (Buchholz et al. 1995). 

Direct comparisons of stability and control derivatives 
identified from flight tests with values identified from 
simulation math models can be used to derive correction 



(b) 

Figure 13. XV-15 tilt-rotor aircraft; (a) hover configuration, 
(b) cruise configuration. 




Figure 14. XV-15 tilt-rotor simulation model validation for 
170 kts. 

factors for significantly improving the model fidelity. For 
example, 
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Figure 15. Tilt-rotor math model error functions and pro- 
posed fidelity criteria for Level D simulators. 

^corrected = f(nonlin.sim.eqns) 

+ [(Lp)flight - (Lp)simJ P 
+ [(Lfi,., > flight " ( L 5 lat )sim ] 5 lat +••• 

( 2 ) 

Identification tools provide a systematic and accurate 
approach to determine these correction factors which are 
routinely used by the simulator industry to improve 
dynamic fidelity. 

Comprehensive simulation studies are often used to define 
flight-control system hardware requirements such as actu- 
ator and sensor filter bandwidths. Franklin et al. (1991) 
used CIFER® to determine actuator bandwidth require- 
ments for a conceptual A/STOVL aircraft. The spectral 
characteristics of the stabilization and command augmen- 
tation system (SCAS) commands to the aircraft control 
surfaces were obtained for ensemble analyses of simulated 
flight tasks, and are shown in figure 16. The results indi- 
cate a SCAS command signal bandwidth (frequency at 
-3dB amplitude) of about 4 rad/sec in pitch and roll, with 
a significantly lower command bandwidth for the thrust 
(vertical) axis. Actuator hardware response bandwidths 
should be 5-10 times the respective SCAS command 
bandwidths to avoid introducing significant phase lag in 
the control loops (Franklin et al. 1991). Similar analysis 
techniques were used by Blanken to determine the effect 
of control-system design on changes in pilot control band- 
width (workload) and handling qualities. This study 
included an interesting comparison of pilot work 



Figure 16. Power spectra of ASTOVL propulsion system 
commands. 

load for the simulation and flight-test environments 
(Blanken and Pausder 1994). 

6. Development 

At the development stage, flight-control system hardware 
and software components and subsystems undergo bench 
testing to verify that the performance characteristics meet 
the design specifications. Sophisticated flight-control 
development facilities (DF) or “hot-benches” allow the 
test of prototype flight software and hardware integrated 
with the simulated aircraft dynamics. In helicopter devel- 
opment, model or full-scale rotors are dynamically tested 
in the wind tunnel and the responses are validated against 
design requirements and comprehensive analysis models. 
This section presents system-identification techniques to 
support development stage validation. Examples are 
drawn from the NASA VSRA project, helicopter actuator 
tests, and the Sikorsky Bearingless Main Rotor (SBMR) 
full-scale rotor wind-tunnel tests. 

An extensive development facility has been used in the 
NASA vertical/short takeoff and landing (V/STOL) sys- 
tems research aircraft (VSRA) project, which equipped a 
YAV-8B Harrier aircraft with a fly-by-wire research 
flight-control system (Foster et al. 1987) (fig. 17). The 
overall flight-control goals of the VSRA program are to 
assess critical technology elements for advanced short 
takeoff/vertical landing (STOVL) aircraft, including: inte- 
grated flight/propulsion control, advanced control and 
display laws, and reaction-controlled bleed-flow require- 
ments. The role of the DF has been for verification of 
control-law flight software, system software, and safety 
monitoring. Actual flight computers and flight hardware 
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Figure 17. NASA V/STOL system research aircraft 
(VSRA). 


were included in the DF to validate flight systems during 
the final development stage. The aircraft dynamics are 
simulated by the VSRA math model, with inputs from a 
test console or a rudimentary pilot-cockpit station. 
CIFER® was exercised extensively to validate broken- 
loop and end-to-end closed-loop frequency responses of 
the DF flight systems against the design models and theo- 
retical analyses. Signal processing and conditioning algo- 
rithms and digital timing were also verified during DF 
testing. 

Actuator system dynamics comprise an important compo- 
nent of the overall high-frequency phase lag in modem 
flight-control systems. Therefore, flight-control system 
stability margins and overall closed-loop performance and 
handling qualities can be significantly degraded if the 
actuator dynamics do not meet the design specifications. 
System-identification bench testing of aircraft actuators 
ensures that expected performance is achieved and that 
costly modifications can be avoided at the flight-test 
stage. Frequency-response identification and transfer- 
function modeling from a typical helicopter actuator test 
are shown in figure 18. Excellent coherence is achieved 
over a broad frequency range (0.2-40 rad/sec) using a 
computer-generated frequency-sweep excitation. The 
actuator dynamics are well characterized by the damped 
second-order response obtained by CIFER® (fig. 18). 
These component system results are used to update simu- 
lation math models and to optimize flight-control system 
gains prior to first flight. 

Structural analysis programs such as NASA structural 
analysis (NASTRAN) are rarely able to accurately predict 
the flexible response beyond the first elastic modes of a 
new aircraft. Therefore, structurally-scaled models or full- 
scale structural test vehicles are evaluated in special rigs 
to verify the elastic characteristics and make final adjust- 
ments to the structural compensation (e.g., notch filters) in 
the control system prior to first flight. Automated 
test/analysis facilities excite the individual structural 
modes of the aircraft with shakers and then use system- 
identification methods to determine model characteristics. 
For modem rotorcraft development, system identification 
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Figure 18. Helicopter actuator response identification and 
modeling. 


has also been effective in extracting dynamic response 
model of subscaled or full-scale rotor systems from 
dynamic wind-tunnel test data. The control response 
dynamics of the SBMR were determined in a joint 
NASA/Sikorsky test in the Ames 40- by 80-Foot Subsonic 
Wind Tunnel (fig. 19) (Tischler et al. 1994). Computer- 
generated frequency-sweep excitation signals to the 
SBMR swashplate actuator were carefully designed to 
ensure adequate identification within the limitations of the 
rotor and wind-tunnel stand. Rotor blade and hub moment 
frequency responses were then extracted using CIFER® 
and were compared to comprehensive simulation models 
of the SBMR. CIFER® was also used to extract the 
rotor’s physical parameters based on a linearized 14 DOF 
analytical formulation of the SBMR dynamics (Tischler 
etal. 1994). 

Figure 20 shows the identified on-axis roll moment 
response to a lateral stick input. The simulation math 
model and 14 DOF identified model agree closely with 
measured responses. The off-axis pitching moment to 
lateral stick input is shown in figure 21. Here, the simula- 
tion model phase response deviates significantly from 
both the measured response and identified parametric 
model, indicating a poor prediction of rotor cross- 
coupling. 
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The key identified physical parameters of the rotor system 
are compared with the GenHel simulation values in 
table 3 (both are updated from the earlier results of 
Tischler et al. 1994). Many of the important rotor parame- 
ters such as Lock number, blade inertia, and effective 
hinge off-set compare very favorably, and reflect the good 
on-axis response prediction of the simulation model. The 


important difference between the identified model and the 
GenHel simulation is the control-phase angle. This 
parameter has a known geometric value of -14 deg in the 
wind-tunnel tests, but the identified value needed to cap- 
ture the measured off-axis response of figure 21 is 
-23.4 deg. This discrepancy in control phasing indicates a 
fundamental problem in the aerodynamic modeling of the 
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rotor. Follow-on analysis of these results have yielded a 
new approach for correcting the simulation math model 
and improvements in the identification methods for free- 
flight results (Takahashi et al. 1995). Accurate cross- 
coupling prediction is especially important for the design 
of decoupling compensators in modem rotorcraft flight- 
control systems. Corrections to the flight-control laws 
prior to final flight software installation and vehicle test- 
ing reduces development flight-test costs and improves 
the final performance of the system. 

7. Flight Testing 

The flight-test program for flight-control and handling- 
qualities validation and optimization has a significant 
impact on the overall development schedule and cost for 
modem fly-by-wire aircraft. System identification pro- 
vides a critical technology for tracking aircraft dynamic 
response performance into flight, solving problems that 
arise in flight tests, and rapidly optimizing control system 
parameters. This section presents system-identification 
methods for control system flight testing. Right data 
results are presented for the VSRA and UH-60A Rotor- 
craft Aircraft Systems Concept Airborne Laboratory 
(RASCAL) projects. 

Right test verification of aero-servo-elastic stability 
margins is an important concern for modem fly-by- wire 
aircraft, where dynamic coupling of the high-gain flight 
control system with light-weight structural dynamics can 
degrade flutter stability. Flutter margin verification using 
system identification has been adopted by British 
Aerospace in the development of a series of fly-by- wire 
high-performance aircraft, as described by Caldwell 
(1994). Near real-time system identification was 
employed during the X-29 aircraft flight testing (Clarke 
et al. 1994) for on-line verification of stability margins in 
a highly-efficient flight envelope expansion program. 
Piloted frequency sweeps were used to excite the vehicle 
structural modes at each test condition, and the tele- 
metered data were then analyzed using high-speed array 
processing computers. Once the stability margins were 
verified, the pilot was cleared to proceed to the next flight 
condition, avoiding the normally time-consuming test 
technique of clearing one flutter test point per flight. In a 
similar application of near real-time identification 
technqiues, CIFER® was used to support flight tests of 
the “Pathfinder,” a large high-altitude solar Unmanned 
Air Vehicle (UAV) (Domheim 1995). Servoloop stability 
margins were extracted based on telemetered data from 
computer-generated frequency-sweep tests, and then 
compared with simulation predictions. When the CIFER® 
results indicated a loss of stability margins at a high 


altitude flight condition, the ground station pilot executed 
real-time switching commands to adjust the Pathfinder 
control law gains. 

Theoretical analyses of the XV- 15 tilt-rotor aircraft 
(fig. 13) predicted that the reduction of whirl mode flutter 
stability margins with increasing flight speed would limit 
the aircraft’s usable flight envelope. An extensive flight- 
test program was conducted to verify the expected mar- 
gins. Early testing using the traditional dwell-delay 
method proved time-consuming and resulted in consider- 
able data scatter. Acree and Tischler (1993) conducted 
automated frequency-sweep tests using wing flaperon 
excitation and subsequently analyzed the data using the 
CIFER® identification tools. The frequency-domain test 
technique proved to be much more time-efficient, and the 
results showed both a reduction in the scatter at specific 
conditions and an improvement in consistency across 
flight conditions. 

Automated frequency-sweep flight testing was also con- 
ducted on the VSRA YAV-8B aircraft (Foster et al. 1987) 
(fig. 17) to determine the locations of the (open-loop) first 
and second structural wing-bending modes, and to verify 
actuator and sensor processing dynamics. The parametric 
model shown in figure 22 was obtained from CIFER®, 
and includes the rigid-body response and second-order 
representations of the two structural modes. Notch filters, 
included to avoid coupling of the flight-control and aeroe- 
lastic dynamics, and control-law gains were subsequently 
updated based on these identification results. Piloted 
frequency-sweep flight testing was also conducted in the 
VSRA program to document the final stability margins 
and closed-loop response for a number of flight condi- 
tions. The broken-loop pitch response for 120 knots as 
obtained from CIFER® is shown in figure 23. The figure 
shows that the dynamics are conditionally stable, with a 
minimum crossover frequency of 1 rad/sec required for 
closed-loop vehicle stability. The nominal crossover fre- 
quency of 4 rad/sec yields a phase margin of 40 deg 
(acceptable). A gain margin of about 8 dB is indicated 
over the broad frequency range 15-30 rad/sec were the 
phase curve has a nearly constant value of about 
-180 deg. 

The identified closed-loop response dynamics are shown 
in figure 24. In the frequency range of 0.3-5 rad/sec the 
response is accurately modeled by a well-damped second- 
order system: 


e 6.35 e -0 - 0485 * 
e^“ [0.953,2.47] 
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Figure 22. YAV-8B aero-elastic wing-bending identification 
(120 kts). 


These dynamics closely match the design response of: 


e com [i.o, 2.0] 

The small equivalent time delay of x = 48 msec reflects 
the VSRA high-bandwidth fly-by-wire actuators and rapid 
digital calculations, and suggests no time-delay related 
handling-qualities problems. 

In some applications, simulation math models are not suf- 
ficiently accurate for control-law design prior to first 
flight. For example, the current state-of-the-art of rotor- 
craft flight dynamics simulation yields a fair prediction of 
the on-axis characteristics, but usually an inadequate pre- 
diction of the cross-coupling response as in the baseline 
simulation result of figure 21, often not even correct in 
sign (Curtiss 1992). Rotorcraft math models are still use- 
ful for initial simulation and control-law development 
efforts, but are less satisfactory for the final determination 



Figure 23. VSRA broken-loop pitch response (120 kts). 


Flight data 



Figure 24. VSRA closed-loop pitch response (120 kts). 


of stability margins and decoupling controller gains. 
Initial flight tests with the SCAS-OFF or with reduced 
control system gains can be conducted to identify new 
aircraft dynamics or to update the simulation for final 
control-law parameter selection. The DLR developed a 
high-bandwidth flight-control system for the Bo- 105 
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variable-stability aircraft (ATTHeS) based directly on 
b are-airframe state-space models extracted from flight 
data using frequency-domain system identification (von 
Grunhagen et al. 1994). This direct use of flight-identified 
state-space models for control-law design represents the 
most sophisticated and demanding application of system- 
identification tools. 

An approach similar to that of the DLR has been adopted 
by the AFDD/NASA in the development of an advanced 
fly-by-wire flight-control system for the RASCAL 
UH-60A helicopter (Takahashi et al. 1995), which uses 
the same airframe as the ADOCS demonstrator (fig. 3). 
Extensive theoretical studies of combat rotorcraft control- 
law concepts for application to RASCAL have been con- 
ducted by Takahaski (1994) and Cheng et al. (1995) based 
on UH-60A simulation math models. At the same time, 
Fletcher (1995) has conducted UH-60A flight tests and 
comprehensive frequency-domain identification studies to 
extract high-order state-space models of the aircraft for 
hover and cruise flight conditions. These efforts were 
brought together in the RASCAL control-law study 
described in Takahashi et al. (1995). Figures 25 and 26 
compare two flight-mechanics-simulation math models 
(“A” and “B”) used for the control-law designs with the 
bare-airframe flight-test data. The on-axis roll response 
agreement between the math models and the flight-test 
data is reasonable at mid-frequency (0.8-10 rad/sec), but 
is inadequate beyond 10 rad/sec due to errors in the pre- 
diction of the in-plane rotor response. Large errors are 
also seen at low frequency. The simulation models show 
poor predictive capability for the cross-coupling response 
of roll rate to longitudinal stick input, with large phase 
errors in the critical frequency range of 1-10 rad/sec 
(fig. 26). While the simulation models were sufficient for 
the preliminary flight-control and simulation studies, they 
are clearly inadequate for selecting final flight gains — 
especially for the response decoupling parameters. 

The identified higher-order linear model is compared with 
the flight data and the simulation models in Figures 25 
and 26. Significant improvement in the on-axis prediction 
is seen for both the high-frequency (rotor response) and 
lower-frequency dynamics. The identified model also 
tracks the off-axis magnitude and phase very closely, 
showing clear improvement compared to the two simula- 
tion models. The excellent predictive capability of the 
identified model is also seen in the time response compar- 
ison of figure 27. 

The identified state-space model was then substituted into 
the model-following control system block diagram in 
place of the original simulation model (“A”) response to 
check the expected flight characteristics. Figure 28 shows 
The identified state-space model was then substituted into 




Figure 25. UH-60A on-axis roll-rate response to lateral 
stick (hover); comparison of simulation and identified state- 
space model with flight data. 



Figure 26. UH-60A off-axis roll-rate response to longitu- 
dinal stick (hover); comparison of simulation and identified 
state-space model with flight data. 
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Figure 27. Time response comparison of UH-60A identi- 
fication model and flight data. 
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Figure 28. RASCAL UH-60A broken-loop roll response in 
hover. 


the model-following control system block diagram in 
place of the original simulation model (“A”) response to 
check the expected flight characteristics. Figure 28 shows 
that the design phase margin is significantly degraded 
when the identified model is incorporated. Further, the 
level of closed-loop cross-coupling (fig. 29) increases by 
20 dB (a factor of 10) in the critical handling-qualities 
frequency range of 1-10 rad/sec. The control-system 
design parameters were then retuned for the identified 
model response. Figure 28 shows that the original design 
crossover frequency, phase margin, and gain margin are 
recovered. Also, the cross-coupling level for the retuned 
system closely tracks the coupling levels for the original 
control-system design (fig. 29). 


ID model 



Figure 29. RASCAL UH-60A coupling response of roll-rate 
due to longitudinal input for hover. 


The full exploitation of system-identification tools early 
in the flight-test development and control system opti- 
mization effort has been illustrated for the Bo- 105 
(ATTHeS) and UH-60A (RASCAL) programs. This 
approach will significantly reduce flight-test development 
time for new aircraft, and will expedite the optimization 
of flight-control system performance and handling 
qualities. 


8. Concluding Remarks 

1. System identification is a full life-cycle technology that 
supports aircraft flight-control system development from 
design specification through flight-test optimization. Sig- 
nificant reductions in development time and costs are real- 
ized by tracking open and closed-loop dynamic response 
characteristics through the development process. 

2. Frequency-domain system-identification methods are 
well suited to aircraft flight-control development since 
many current design specifications, design and analysis 
techniques, and acceptance flight-test techniques are 
based in the frequency domain. 


17 


3. Reliable computational tools for system identification 
are available and have been successfully employed in 
many recent aircraft programs. 

4. System identification is especially effective in provid- 
ing a transparent and integrated understanding of 
handling-qualities characteristics and system stability. 
Considerable improvements in system performance are 
facilitated by the rapid availability of accurate end-to-end 
and subsystem dynamic models. 
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